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Abstract 

Metabolic allometry, a common pattern in nature, is a close-to-3/4-power scaling law be- 
tween metabolic rate and body mass in organisms, across and within species. An analogous 
relationship between metabolic rate and water volume in river networks has also been observed. 
Optimal Channel Networks (OCNs), at local optima, accurately model many scaling properties 
of river systems, including metabolic allometry. OCNs are embedded in two-dimensional space; 
this work extends the model to three dimensions. In this paper we compare characteristics of 
3d OCNs with 2d OCNs and with organic metabolic networks, studying the scaling behaviors of 
area, length, volume, and energy. In addition, we take a preliminary look at comparing Steiner 
trees with OCNs. We find that the three-dimensional OCN has predictable characteristics anal- 
ogous to those of the two-dimensional version, as well as scaling properties similar to metabolic 
networks in biological organisms. 

1 Introduction 

Fractal Networks have been widely used to study well known network behaviors. Some of the 
properties commonly used in the study of fractal networks are scaling factors (or laws); scaling 
laws can describe how lengths are related to area and areas are related to volumes. In this paper, 
we further study these allometric scaling laws. We use simulation and elementary mathematical 
bounding techniques to study such laws. 

A common pattern in nature, known as metabolic allometry, is a 3/4-power scaling law between 
metabolic rate (Af) and body mass (B) in organisms, across and within species: B cx An 
analogous relationship between water mass (C) and contributing area (A) in river networks has also 
been observed: A cx C^^^. 

The relationship in biological organisms was first observed by Max Kleiber in 1932 and has 
garnered a lot of attention in the last decade and a half from biologists, mathematicians, and 
physicists alike |13| . The areas this relationship finds application in include models of climate change 
effects in rivers |19j and oceans [l5j . analyses of river stability [17 , and in studies of information 
networks |18j . Many theories concerning the origins and nature of this phenomenon have been 
proposed, and there has been much debate about how ubiquitous the scaling law really is. 

In 1997, West, Brown and Enquist (WBE hereafter) introduced a theory purporting to explain 
metabolic scaling [2^1 . Their theory is based on the idea of a fractal space-filling hierarchical network 
of resource distribution, with certain assumptions based on simplified biology. The geometry of the 
network (e.g. branching ratios, vessel radius relationships) are important to the development of 
their result. The authors claimed the model is applicable even to organisms without a physical 
blood vessel network since they can be treated as having a virtual distribution network |27| . 

Subsequent studies questioned the WBE results, by showing that data does not appear to strictly 
follow a 3/4-power scaling law after all; the theoretical approach taken by WBE has also been 
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criticized [IlllH]. Some authors claim that an exponent of 2/3 fits the data better [5], while others 
assert that no universal scaling law exists at all [TJ- Savage et al. pointed out that the WBE model 
predicts a curvilinear relationship between logB and logM, and only predicts a constant power-law 
relationship in the infinite limit of body size; however, even with finite size adjustments and variations 
in the structural properties of the model, they found the results were still inconsistent with data 
pS] , Some other modifications to the WBE model do yield results somewhat more consistent with 
data, such as taking into account body temperature [S] or fluid velocity [1]. Kolokotrones et al. use 
regression on large set of organism data to show a quadratic, rather than linear, relationship between 
logB and logM [M]. Also, their results reported a lot of variation. They build off the modifications 
to WBE in [5S] and find that a model network with a proportional transition between area-increasing 
and area-preserving branching yields a fit to the data almost as good as their empirically-derived 
quadratic model. 

Other theories and models have been put forth, some based on network properties like the WBE 
model, and others on other biological properties. Dodds [10] proposed a model based on virtual 
networks of metabolite transportation rather than physical networks. 

Banavar et al. proposed a general model for any efficient transportation network, predicting 
B (X M^/(^+^^ (at most) where D is the dimension of the network (i.e. D = 3 for organism 
metabolism, and D ~ 2 for river networks) [3]. Thus, they predict that if the network were as 
efficient as possible, organism metabolism should scale as i? oc M^/'*, and deviations of the exponent 
below 3/4 might be explained by some inefficiency in the metabolic distribution network. Later the 
authors added to this model a supply and demand mechanism that accounts for the deviations 

Isaac et al. propose a less strict understanding of the "3/4-power law": though the 3/4 exponent 
might not be a consistent, universal constant, it is still an observable trend; if there isn't a single 
unifying principle explaining why it appears, there still may be multiple reasons for it worth studying 

m- 

In the context of river networks, the analogous properties to metabolic rate and body mass are 
contributing area (A) and water mass (C). It is well-observed that the properties obey a power-law 
relationship with an exponent of slightly more than 2/3: A oc C^/'^. This is consistent with the 
general resource transportation model of Banavar et al. ^ for systems in two dimensions, and also 
with the virtual network model of Dodds 10 . 

Rinaldo et al. devised a model of river networks called Optimal Channel Networks (OCNs) pi] . 
This model uses a grid network to represent the area of a river basin, and constructs a spanning 
network on the grid that minimizes a functional representing the energy required to drain the area 
of land through river channels. The energy functional is derived both from theoretical physical 
characteristics and observed characteristics of river basins [23]. Statistical properties of Optimal 
Channel Networks (including metabolic allometry) correspond excellently with observed properties 
of real river networks [151 [55] . 

This paper explores a three-dimensional version of OCNs, using simulations to compare char- 
acteristics of 3d OCNs with 2d OCNs and with organic metabolic networks and study the scaling 
behavior of area, length, volume, and energy. We do not attempt to explain 3/4-power scaling in 
nature, but rather explore some interesting and potentially useful connections, somewhat in the 
spirit of Isaac et al.'s proposed understanding of the scaling law [T5] . 

We further study an alternative model using Steiner trees for OCNs and we compare the scaling 
behaviors between these two models on small network sizes. 

Section 2 describes definitions and methodologies used in our paper. In Section 3, we derive 
lower bound results for various scaling laws. In section 4, we describe our simulation results and 
show Energy, Length and Volume Scaling with respect to area and compare it with the lower bound 
results described in Section 3. In Section 5, we describe an alternative model based on Steiner Trees 
and derive analytical results for small network sizes. Section 6 gives conclusions and suggestions for 
future work. 
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2 Definitions and Methods 



Consider an area of land, and embed a grid network in it such that every node in the network is 
associated with a unit area of land. Construct a spanning tree on the grid and direct it so that all 
nodes have a path in the tree to a single root (call it the outlet). Define a partial order on the nodes 
such that a; < y iff a; is on the path from y to the root. If a; < we will say y is upstream of x. The 
terms area, volume, upstream length, link length, and energy are defined here in this context. 

Note that each node in the spanning tree, except for the root, has exactly one link directed out 
of it. Thus we can name each link according to its corresponding node; that is, when we speak of 
link X we are referring to the link out of node x. All the properties defined here for nodes will be 
extended to their corresponding links. For example, the area of link x is defined to be the area of 
node X. 

The area of node x, A^, is the total number of nodes whose paths to the outlet include x] that 
is, the number of nodes y such that x < y. This is equivalent to the sum of the areas of each node 
directly linked to x plus one for x itself: 

A, - Y.Ay + 1 (1) 

y 

where y ranges over the nodes directly linked to x [2T] . 

The volume of node x, Cx, is the sum of all the areas of all nodes upstream of a; [3]: 

— ^ Ay (2) 

y\x<y 

Equation |2] can be transformed into a recursive definition for volume: 

y 

where y ranges over the nodes directly linked to x. 

Upstream length is another quantity of interest; it is defined as the number of links in the path 
from a node a; to a source node, where at each step upstream the path taken goes to the node with 
the largest area [71. Statistically, it is equivalent to define it as the distance (counted in number of 
links) to the farthest away upstream source [20] . 

Link length is a weight assigned to each link, generally representing the geometric length of the 
link in the embedding. The length of link i is denoted by li. 

Energy is a function of the link lengths and areas of all the links. The energy of a given tree 
configuration s is given by: 

H,{s)= ^^'^ (4) 

i^links{s) 

with < 7 < 1- A tree network that achieves the minimum Hj(s) over all possible tree configurations 
on a given grid is an optimal channel network |21| . 

Figure [l] shows a grid and a possible tree configuration si. si minimizes iJi/2(s) for all config- 
urations on this grid. Area, volume and upstream length values are also depicted for the nodes in 

Si. 
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(a) Underlying grid 



(b) Areas of nodes 
in si 



(c) Volumes of nodes 
in si 



(d) Upstream lengths 
of links in si 
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Figure 1: The outlets are in the top left corners. Tree si minimizes Hi/2{s) for the grid in 
-H"i/2(si) = 24.3024964149. 



Rinaldo et al. [H] derived a gamma value of 1/2 from estimations of physical properties of rivers, 
and hypothesized that natural rivers would minimize Hi/2{s). A remarkable result of |21j is that 
while the global optima of OCNs do not match well statistically with real river networks, the local 
optima are an excellent match. The idea that nature tends towards a local optimum makes sense, 
so this suggests that OCNs are a good model for what actually happens in a river [55] . 

When 7 = 1 and all links are of unit length, the energy is equivalent to the volume. When 7 = 0, 
the energy is equivalent to the total length of all links, so minimizing the energy is equivalent to 
finding a minimum spanning tree (which, when all links are of unit length, would be any spanning 
tree). 

In most studies of OCNs, the network of potential links is a 2-dimensional lattice with additional 



links between diagonally-adjacent pairs of nodes (Fig. 2a). For the sake of simplicity, all the links, 
diagonal and orthogonal, are here considered to have unit length, because it has been shown ([23 ) 
that the properties with which this paper is concerned are independent of whether the diagonal links 
are given realistic length or not. 





(a) 7x7 network of potential links 



(b) 2x2x2 3d network of potential links 



Figure 2 

In this study we also analyzed a three dimensional version of optimal channel networks, built on 



a 3-dimensional grid (Fig. 2b I . 

To make an OCN, we first generated a random spanning tree on a given grid. We used Prim's 
algorithm for minimum spanning trees for this. Essentially, the algorithm builds the tree from a 
root node, maintaining a set of frontier edges from the grid that are incident on exactly one node in 
the tree. Our algorithm randomly chooses one edge at a time from this set, adds it to the tree, and 
updates the set. Rinaldo et al. also began with random spanning trees, though they did not specify 
their method for generating them |21j . 

The random tree was then optimized using a version of Lin's algorithm for TSP, as in [21]. That 
algorithm is as follows: 
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Input: Undirected graph G; random spanning tree T, directed so that there exists a path in 

T from each node to a single root 
Output: Locally optimal spanning tree on G 
while Convergence has not been reached do 
Step 1: Choose a random node i in T 

Step 2: Choose a random neighbor of i in G distinct from the node i 
links to in T 

Step 3: Construct T* from T by redirecting i to link to this new node 
if T* contains a loop then discard T* 

else if H-^{T*) < H^{T) then discard T and let T* be the new T 
else discard T* 

The condition for convergence was that the number of improvements made be only 1% of the 
total number of iterations (or 2% for larger networks). See Figure|3]for an illustration of the process. 




(a) Initial tree configuration T. The outlet is in the 
top left corner. H^{T) = 9.65028. 




(c) Step 3: A new link for i in the tree has been chosen, 
creating T* . H-,{T*) = 9.8637 > H-,(T), so this tree 
will be discarded. 




(b) Step 1: Node i has been chosen; the unused links 
to its neighbours in G are highlighted. 




(d) the next iteration, Step 3: A different new link for 
i in the tree has been chosen, creating T* . H^{T*) = 
9.44949 < Hj{T), so this tree will become the new T. 



Figure 3: Optimization algorithm 
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Figure 5: An OCN on a 3x3x3 grid; the red hnks arc in the OCN. 

Figure [4a| shows a random tree on a 60x60 grid and Figure [4b| shows the result of using the above 
algorithm to optimize the tree. In order to make the structure of the network clearer, links are 
drawn with a thickness proportional to the log of their area. Figure [5] depicts a three-dimensional 
OCN on a 3x3x3 grid. 

A range of grid sizes were analyzed. Two-dimensional rtxri grids with n ~ 10, 20, 30, 40, 50, 60, 
70, and 80, and three-dimensional nxnxn grids with n = 8, 10, 12, 14, 16, 18, and 20 were used. 
For the smaller grids, data was averaged from at least 4 realizations of each size, while for the larger 
grids at least 3 realizations were used for each. The quantities measured were distributions of area 
and length in whole basins, and average length, volume, and energy per area of subbasins and whole 
basins. 

3 Analytical Proofs 
3.1 Energy bounds 

In [S], Colaiori et al. proved lower bounds for the energy of OCNs on an orthogonal grid (one with 
no diagonal links). Here, we look at lower bounds for the energy of OCNs on an eight-neighbor grid. 
The main result for this section (Theorem 3) is that for an n x n 8-way grid, — |n -I- 1 is a lower 
bound on Hq^{s). 

Lemma 1. Let (P) he the optimization problem 

n 

min AJ 

i=l 
n 

s.t. Ai = n + m 

i=l 

Ai>iyi 
A,€Zyi 

where n is the number of entries in the vector A, m is some natural number, and < 7 < 1. 

Al = A2 = ... = An-i = 1, An = m + 1 is an optimal solution to (P). 

Proof. Suppose A is an optimal solution for (P) and 3j s. t. 1 < Aj < m + 1. 
3k s.t. 1 < Ak < m + l,j ^ k. 

We can construct A* such that A* = Ai Vi 7^ j, fc; A* = 1; = Aj + Ak-l. 
Then Er=i AV = (EIli ^7) -A]-Al + l + {A, + A^ - 1)^. 

Since f{x) = when < 7 < 1 is concave, and Aj^A^ > 1, {Aj + Ak - 1)'^ < A'J + Al - I. 
So J27=iAV — E"=i^7- Thus A* must be an optimal solution for (P). By an iteration of this 
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argument, there must be an optimal solution with all entries except for one equal to 1. Without 
loss of generality, the entries can be rearranged so that is the largest (only non-unit) entry. □ 



Theorem 2. For an n x n eight-neighbor grid G, the optimal -ffi(s) is 



C\ C'2 (.3 C'4 (.5 

Figure 6: The five "stripes" of a 5 x 5 8-way grid with the outlet in the bottom left corner. 

Proof. Consider the subsets of vertices Ck where vertex i g iff the distance from i to the outlet 
equals fc — 1 (Figure |6]). In an n x n 8-way grid there are n such stripes. Assuming the outlet is 
located in a corner, \Ck\ = 2k — 1. Clearly, 



E ^7 = EE^7 



(5) 



k=2 ieCk 



Recall that the area of node i is the number of nodes upstream of i, plus 1 for i itself. So for a 
given fc, X^ieCfc ^7 ^^^^ have to include the number of nodes in Ck and the total number of nodes 
upstream of all i G Ck, which must include all the nodes j G U;>fe C';, since their paths to the outlet 
must pass through stripe Ck at some point. C; is simply the set of all the nodes in the grid 

minus those within the fcxfc section enclosed by Ck, which is — fc^ nodes. Thus the total amount 
of area being passed through Ck is at least 



\Ck\ 



(6) 



Then we have 



^i(^) - EE^' ^ E(2fc-i 



fe=2 ieCfc 



k=2 



The rightmost sum simplifies to |(4n'^ — in? — n); therefore this is a lower bound on Hi{s). 

Consider the network t where every link from node i to node j is such that ifi € Ck then j G Ck-i- 
Here, quantity pi is exactly the area passing through stripe Ck, so Hi{t) = |(4n'^ — — n). Hence 
the lower bouna is achievable. 

□ 



Theorem 3. For an n x n eight-neighbor grid, — ^ti -f- 1 is a lower bound on Hq c^{s 
Proof. Following the same reasoning as above, we have that 

n 

^7(-) = EE^z 

k=2 ieCk 



(7) 



7 



and the total amount of area being passed through Ck is at least |Cfc| + — fc^. 

By Lemma [l] the optimal way to distribute this area over the vertices in Ck is to send all the 
upstream area, — k"^ , through one vertex in Ck, leaving each of the other vertices in Ck with area 
1. (Note that this is not usually feasible, but it is a lower bound.) 
Thus for a given k, 

Y,A'l>\Ck\-l + {n^-k^ + lT 

Since f{x) = is concave, 

n n n 

k=2 k=2 k=2 

So 

n 71 n 

^ ^ > Y.[2k -2 + {n^-k^ + ly] > n^-3n + n^^+' -n^^~Yl 

fc=2ieCfc fc=2 k=2 

For 7 = 0.5, this evaluates to ^ri^ - An - - 1) = f "-^ - + 1. □ 



3.2 Volume Scaling 

In |16| . Maritan et al. make an analytical prediction for the scaling of volume with area in rivers. 
Their analysis depends on the relationship {Lx) (x A'^, where (L^) is the mean distance from vertices 
upstream of x to x, and h is Hack's exponent, from Hack's law (a power-law relationship between 
basin length and area) This is combined with the equation Vx — Ax{Lx) to arrive at the 

relationship oc Al^^. This relationship is verified by their data on real river networks as well as 
OCNs, where (for both of which) h is typically close to 0.57 |16j . 

If h were to equal exactly 0.5, this would imply isometric scaling of length and area in river 
basins, as A^ would scale directly with l^, and shape would be preserved. This would line up with 
the prediction in [3j that V cx Ax^^^^ in the most efficient networks for D — 2. However, since 
river basins tend to elongate with growth, getting proportionally narrower as they get larger, h is 
usually slightly greater than 0.5, and the lower bound predicted in [5] is rarely reached. 

Shifting focus to three dimensions, we can look for evidence of a corresponding Hack's law for 
three dimensional OCNs (that is, evidence that area scales with length to some constant power). 
If there is a relatively constant scaling exponent h, the analysis in |16j can be extended to three 
dimensions, since the steps taken there do not otherwise depend on dimension. Thus, we might 
expect to see Vx oc A].'^^ for some h; according to the analysis in [3], 14 should scale at least as 
At^^ , so we might expect 1/3 to be a lower bound for h, with h slightly greater than 1/3 if three 
dimensional OCNs elongate in a similar way to two dimensional ones. These predictions are borne 
out in the simulation results that follow. 
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4 Experimental Results 
4.1 Energy Scaling 
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Lower Bound 
-Linear fit: f(x) = 



0.786x'^1.158 R^ = 0.9991904447 



♦ Observed Values 
— Power law fit: 



f(x)= 1.482X- 64.061 
R== 0.9999795166 




10000 



100000 



Area 



Figure 7: Log-log plot of area vs. energy for 2-dimensional OCNs, showing both analytical lower 
bound (lower, blue) and observed values (upper, red). 

Fig. [T] shows how empirically observed minimum possible energy scales with area in two-dimensional 
OCNs. The estimated line of best fit is shown, as well as the analytical lower bound derived in Section 
[Sj The observed values of energy follow a power law with a small exponent, while the lower bound 
is essentially linear with respect to area. 



4.2 Length Scaling 

Scaling exponents h where / cx in two and three dimensions is shown in Table [ij The observed 
values of the exponents for two dimensions are within the bounds found in other studies [111 116) . 
For three dimensions, the exponent is also fairly consistent, though less so than in two dimensions. 
It also deviates slightly more from the isometric value of 1/3 than the two dimensional version does 
from 1/2. The fact that it is higher than 1/3 fulfils the hypothesis that three-dimensional basins 
elongate like two-dimensional ones. 



basins analyzed 



2d whole basins 
2d all subbasins 
3d whole basins 
3d all subbasins 



0.6066895125 
0.5808520176 
0.3645961454 
0.4230582779 



Table 1: Length scaling exponents 



4.3 Volume Scaling 

Fig. [8] shows the scaling of volume with area in two-dimensional OCNs, for whole networks and for 
all subbasins. Note that the exponents are very similar, indicating that the scaling behavior is the 
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same within basins as across different sizes of basins. The exponents are both close to 1.57, as in 



600000 ■ 
500000 ■ 
400000 ■ 
300000 ■ 
200000 ■ 
100000- 
0- 



(X) = 04936201486 x^l 5701068248 




R== 0.9999109801 









1000 2000 3000 4000 6000 6000 7000 

A 




(a) 



(b) 



Figure 8: Volume plotted with area for whole nxn basins ( 8a ) and all subbasins ( 8b ) of 2-dimensional 
OCNs. 
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1000 2000 3000 4000 5000 6000 7000 8000 9000 
A 



(a) Whole nxn basins 



(b) all subbasins 




f(x) = 04077052272 x^l .4445624572 
10000 R!= 9951362S6 




(c) subbasins of 10x10x10 OCNs (d) all subbasins with Area < 1000 

Figure 9: Volume plotted with area for 3-dimensional OCNs. 
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basins analyzed 












all subbasins 




1.424 


all subbasins of 20x20x20 network 




1.437 


subbasins with area > 1000 




1.3776 


subbasins of 20x20x20 network with area 


> 1000 


1.342 


all subbasins with area < 1000 




1.44456 


subbasins of 20x20x20 network with area 


< 100 


1.457 


subbasins of 10x10x10 network 




1.4506 



Table 2: Volume scaling exponents for three-dimensional OCNs 



The three-dimensional results for volume scaling are displayed in Fig. [9]and summarized in Table 
[2j The differences in exponents between the different data sets show higher scaling behavior in the 
sets of smaller basins; the difference between the exponent for whole basins (9a I and all subbasins 
( |9b[ ) is due to the greater concentration of smaller basins in the set of all subbasins, skewing the slope. 
Note that the exponents in 9c and 9d are similar, showing that scaling behavior in the subbasins is 
the same for the same distributions of area independent of the size of the encompassing whole basin. 
In Section |3] we predicted that a would equal 1 -I- h, using the analysis in |16j . and this relationship 
is clear in Tables [T] and [21 

The exponent a ranges from «1.34 to «1.46, with a higher exponent for collections of smaller 
basins. These results are not inconsistent with the data for metabolic allometry in biological or- 
ganisms, where a power-law fit to data from a wide range of sizes approximates that metabolic rate 
scales with mass to the power of 0.70 (l/a) [25]. Not only is this value within the range found in the 
OCNs, but in organisms as well it is found that sets of smaller organisms yield a higher exponent 
than sets of larger ones [TH [23] . 



5 Alternative Models 
5.1 Steiner Tree Model 

In addition to studying OCNs, we took a preliminary look at comparing OCNs with Steiner trees. 

Steiner trees are minimum-weight trees that connect a given subset of nodes in a graph. The 
specified nodes are called terminals. Steiner trees may include non-terminal nodes; these are called 
Steiner points (or Steiner nodes). 

In [6], Steiner trees built on a grid are considered. The underlying graph is all points in a 
Euclidean plane, completely connected, and the terminals are the points of a grid. We began 
comparing Steiner trees with OCNs by looking at the relative optimality of OCNs minimizing Hq^{s) 
and normal Steiner trees on grids of the same sizes. 

In order to evaluate the energy of a Steiner tree, the tree needs to be directed; this can be done 
by choosing the outlet to be an arbitrary corner node and directing all the edges towards the outlet. 

For the purposes of this study area was defined as the number of terminal nodes in the subtree 
rooted at a given node. This way, the terminal nodes have exactly one unit of area associated with 
each of them, just as the grid nodes do in the OCN. 

The OCNs here have realistic link lengths, where diagonal links are length \/2, to make them 
comparable to the Steiner trees. 

OCNs s of size 2x2, 3x3, 4x4, and 5x5 were constructed and H^{s) for each were compared with 
H^{t) for Steiner trees t of the same size. The Steiner trees used were found in [5]. Figure 10 depicts 



the pairs of trees that were compared. On the left of each pair is the Steiner tree, with open circles 
for the Steiner nodes, and on the right of each is the corresponding OCN. 
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(a) 2x2 (b) 3x3 (c) 4x4 (d) 5x5 

Figure 10: Steiner trees and OCNs for different nxn grids. 



Since Steiner trees minimize total link length, it was clear that they would always have lowest 
energy for 7 = 0. For larger values of gamma, however, where area becomes a bigger factor than 
link length, we hypothesized that the OCNs would have a lower energy value than the Steiner 
trees. Finding whether this was true, and where the crossover point would be, was the goal of the 
comparisons. 

5.2 Analytic Results 




0.02 0.04 0.06 0.08 0.10 0.2 0.4 O.fi 0.8 1.0 



(c) 4x4 (d) 5x5 

Figure 11: Plots of H^{s) (blue) and H^{t) (red), where s is an OCN and t a Steiner tree of the 
given size. In each plot, the horizontal axis is 7 and the vertical axis is energy. 



Figure [TT] shows plots of the energy for corresponding Steiner trees and OCNs, < 7 < 1. As 
expected, the Steiner trees have lower energy than the OCNs for smaller 7, while for 7 past a certain 
crossover point the OCNs are better. Table 3 lists the crossover values of 7 for each size. 
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size 


crossover 7 


2x2 


0.55433153188825812933 


3x3 


0.63036154802337965422 


4x4 


0.01907039679465522977 


5x5 


0.32264708810531986705 



Table 3: Crossover values of 7 



5.3 Interpretation 

From the sizes studied, it is not possible to detect a trend in crossover points. 

The next step would be to combine the two models, constructing Steiner-type trees that minimize 
Hy{s) for varying 7 and comparing them with real river networks and with normal OCNs. 

6 Conclusion 

We have looked at three dimensional OCNs as they compare with metabolite distribution networks 
in organisms. In the organism context, volume is interpreted as volume of blood (which is directly 
proportional to the mass or volume of an organism) and area is proportional to the number of 
capillaries f53]. In this context, the grid used to create OCNs no longer represents fixed space (as it 
does in the two-dimensional/river context) and the amount of physical area directly associated with 
each node may not be the same over different body masses. This would only have a linear effect on 
the relationship between area and mass, however, so this has no effect on the scaling exponent. 

Though the three-dimensional networks do not actually look like the metabolic networks in 
organisms, we have shown that they have similar area/volume scaling behavior. This adds credence 
to the idea that the observed metabolic scaling is a result of characteristics of the distribution 
network. 

Another way to study the possible similarities between OCNs and metabolic networks would 
be to derive an energy functional for three dimensional OCNs from metabolic networks themselves, 
rather than using the one derived from rivers. This could also include consideration of heterogeneity. 
Heterogeneity in the river context is addressed in and [5]; possibly this work could be extended 
to three dimensions. 
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